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Abstract 

The implementation of a two-equation k-co turbulence 
model into the NPARC flow solver is described. Motivation 
for the selection of this model is given, major code 
modifications are outlined, new inputs to the code are 
described, and results are presented for several validation 
cases: an incompressible flow over a smooth flat plate, a 
subsonic diffuser flow, and a shock-induced separated 
flow. Comparison of results with the k-e model indicate 
that the k-co model predicts simple flows equally well 
whereas, for adverse pressure gradient flows, the k-co 
model outperforms the other turbulence models in NPARC . 
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k-co turbulence model constants 
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k-co turbulence model constants 
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Introduction 

In NASA’s High Speed Research (HSR) program, 
computational fluid dynamics (CFD) is increasingly being 
used to design and evaluate inlet and nozzle configurations 
for High Speed Civil Transport (HSCT) applications. 
Accurately predicting flow separations andmixing between 
primary and secondary flows in highly turbulent mixer- 
ejector nozzles as well as pressure losses and bleed flows 
in supersonic mixed-compression inlets is a major 
challenge. The turbulence model employed is often the 
limiting factor in these types of simulations and, as a 
result, is often blamed for providing poor agreement with 
experimental data. Thus, the search for better models is 
ongoing. 

Turbulence models vary from relatively simple algebraic 
models to one-equation, two-equation, and full Reynolds 
stress models. Moving to the next level of complexity 


generally involves solving additional transport equations. 
Thus, as the model complexity increases, so does the 
effort required to implement the model. Higher level 
turbulence models are also more computationally 
expensive because they require more computer memory 
and cpu time. Two-equation models are commonly used 
because they are considered to provide relatively good 
results for the demand they place on computer resources. 

The most widely used two-equation turbulence model, 
the k-e, solves two transport equations for the turbulent 
kinetic energy k and the rate of dissipation of turbulent 
kinetic energy e. Important groundwork on this model was 
performed by Jones and Launder 1 and Launder and 
Spalding . 2 The Jones and Launder model is often referred 
to as the standard k-e model. 

The k-co model is similar to the k-e model except that the 
second turbulent quantity CO is the specific dissipation rate 
(dissipation rate per unit turbulent kinetic energy). 
Kolmogorov 3 was the first to propose this type of model. 
Since then, others have made modifications to the original 
formulation. Among these are Menter 4 and Wilcox 5 . 
Wilcox 6 presents an interesting historical overview of the 
k-co model. 

Recent efforts at NASA Lewis Research Center to 
model the flow field inside subsonic diffusers have 
demonstrated the inability of current models to accurately 
predict flows with adverse pressure gradients. The k-co 
model has been shown to predict separated flows better 
than the standard k-e model, especially in flows with 
adverse pressure gradients . 7 ' 8 As a result, the k-co model 
of Wilcox 5 was recently installed in both the two- and 
three-dimensional versions of NPARC. Because the k-co 
model is a two-equation model, it does not require additional 
cpu time relative to the k-e model, and since it is similar in 
form to the k-e model, a similar algorithm can be used. 
Furthermore, this model has the benefit of improved 
transition simulation and is capable of more realistic 
rough-wall treatments. This report provides an overview 
of some of the turbulence models already available in 
NPARC, a discussion of the implementation of the k-co 
model in the code, and results obtained using the new 
model for several flow fields. 


NPARC Code 

The NPARC code, previously known as the PARC 
code 9 and originally developed at the Arnold Engineering 
Development Center (AEDC), uses the Beam and 
Warming 10 approximate factorization algorithm to solve 
the Reynolds-averaged Navier-Stokes equations at discrete 
grid points. Generalized boundary conditions allow the 
user to specify any portion of a grid as a boundary. 
Noncontiguous multiple block interfacing can also be 
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used to simplify the grid generation of complex geometries. 
Cooper 1 1 discusses recent features of the code. 

The NPARC code contains three algebraic turbulence 
models, the first of which is the Thomas model. 12 It 
computes turbulent viscosity for wall-bounded and free 
shear layer flows and has been optimized for the latter. The 
second, the Baldwin-Lomax model, 13 only computes 
turbulent viscosity in wall-bounded regions. An algebraic 
RNG (renormalized group theory) model 14 is also 
available. 

NPARC also contains the Bald win-Barth one-equation 
model 15 and a two-equation turbulence model based on 
the Chien low-Reynolds-number k-£ model, 16 which is 
presented below in detail to underscore its similarity to the 
k-co model: 

1*1 -Wt (1) 

9(pk) ^ 9(pu;k) 9 f [ m "| 9k 

at 9x, ~ax i [^ + a k Jax i _ 

+n-pe(l + M2)-2n^- (2) 

9(pe) | a(pUj£) _ 9 ( < p t ' 9e 
9t 3xj 9xj ^ c £ J 9xj 

+ c e f l nf-C t2 f 2 Py-2n^-e‘- < '- 5 >'*) (3) 

where 
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f^ = i.o-e<-°- 01,5 y + > 
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f, =1.0 
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f 2 = 1.0-0.22e _(Re ' /6)2 

(7) 

Re,^ 

1 no 
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and C^= 0.09, C e i = 1.35, C £ 2 = 1.80, ak=1.0, a e =1.3. 
The turbulent Mach number is defined as M t 2 = 2k/a 2 , 
where a is the reference speed of sound, and is used in 
Sarkar’ s compressibility correction. 1 7 This compressibility 
correction has been incorporated into the NPARC code to 
compensate for the apparent increase in turbulent 
dissipation at higher Mach numbers. The e-equation (3) 
has been relatively insensitive to compressibility effects 
and therefore requires no modification. 

The terms on the right-hand side of equations (2) and (3) 
correspond to the diffusion, production, dissipation, and 
near-wall damping terms, respectively. One of the 
difficulties with the k-e model is that there is no natural 
boundary condition for £ near a solid surface. Since the 
turbulent kinetic energy and the turbulent viscosity are 
both zero along a wall, the dissipation rate is set to zero as 
well. It is the near- wall damping terms of equations (2) 
and (3) that allow all the turbulent quantities to be set to 
zero near a viscous surface. Another problem with the 
Chien model is that it requires that the distance from the 
wall be calculated. In complicated geometries, and 
especially in three-dimensional cases, this can be difficult 
to compute. Lang and Shih 18 also point out that near-wall 
damping terms which use y + are not desirable, particularly 
near separation regions. 

The Wilcox k-co model does not have the aforementioned 
shortcomings. The value of (0 near a viscous wall can be 
related to the surface roughness of the wall. This condition 
can be enforced at the boundary and does not require the 
use of near- wall terms or y + . The Wilcox k-co model is 
formulated as follows: 

H t =oc*— (9) 

1 CO 
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below, comparisons are made between the models using 
this relation. However, since P* is afunction of Re t , which 
in turn depends upon to, the conversion of £ to co makes the 
assumption that P* = 0.09. 


Code Modifications and Usage 

Because of the similarity in the k-e and k-co formulations, 
the current implementation makes use of the existing 
Chien k-e algorithm available in NPARC. The algorithm 
for solving the k-£ equations dates back to the work of 
Nichols 19 with recent modifications by Georgiadis, 
Chitsomboon, and Zhu. 20 The most significant changes 
occurred in the calculation of the source terms involved in 
the k- and co-equations. The source terms are those on the 
right-hand side of equations (10) and (11) excluding the 
diffusion terms. 

The k-co model is selected by setting IMUTUR (or 
IMUTR2) equal to 7. Since it is similar in form to the 
k-e model and uses the same solution algorithm, users 
should expect stability and convergence to be comparable. 
Users are advised to use appropriate near- wall grid spacing 
along viscous surfaces. This helps to resolve the boundary 
layer and the large co-gradients which occur near these 
surfaces. 

Compressibility Corrections 


II 

CQ. 

(16) 

Re = — 
1 1 1 

(17) 


(ICO 


A compressibility correction has been added to the 
model equations to enhance predictions at higher Mach 
numbers. The corrections described by Sarkar 17 or 
Wilcox 2 1 may be selected by setting the variables COMPK, 
COMPW, and CMT0 according to the following table: 


and 

c k =2 (Xq = p/3 R s =8 R co =2.7 

a ( 0= 2 a 0 = 1/10 R k = 6 

The compressibility correction is given by 

F(M,) = Max(M 2 - M 2 0 ,0.o) 

where the turbulent Mach number is the same as before, 
and the constants default to = 1.0, = 0.0, and 

M t0 = 0.00. 

By comparing the k-co and the k-e models, a relation 
between co and e can be found. Wilcox 4 5 defines this 
relation as co = e/(P*k). In the results to be presented 


COMPK 

COMPW 

CMTO 

Correction type 

1.0 

0.0 

0.00 

Sarkar (k-eq.) a 

1.0 

1.0 

0.00 

Sarkar (k+m eqs.) 

1.5 

1.5 

0.25 

Wilcox 

0.0 

0.0 

0.00 

None 


“Default. 


These variables correspond to ^co.and Mto, respectively, 

in equations (10) and (1 1) and have been included in the 
TURBIN namelist block. As discussed in reference 2 1 , the 
Sarkar correction only affects the k-equation in the k-s 
model. When used with the k-co model, however, correc- 
tions should be added to both the k- and co-equations. 
Computations for the Sajben 22 diffuser strong-shock case 
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(discussed in the Results section) indicate that adding the 
Sarkar correction to only the k-equation provides the best 
results. The Wilcox correction was of little use since the 
turbulent kinetic energy led to an M t value that was below 
the cutoff turbulent Mach number CMTO = 0.25. As a 
result, the Sarkar correction to the k-equation was chosen 
as the default setting, but the other options have been 
included because of the lack of extensive testing of high- 
speed flows. These options should provide the greatest 
amount of flexibility and limit future code modifications. 

Boundary Conditions 


The boundary value for co along a viscous wall is related 
to the surface roughness of the wall. This relation is given 
by Wilcox 5 : 



at y = 0 


where 


S R 


^ 50 ^ 





<25 


4 >25 


(18) 


(19) 


and is valid for kj£ = u x kR/v values up to 400. The NPARC 
code has been modified to allow the user to input the wall 
surface roughness for no-slip walls. This value is read in 
through the auxiliary pressure variable. The variable 
SRDEF has been added to the TURBIN namelist block 
and is used to specify the default smooth-wall surface 
roughness. When this variable is used, the surface rough- 
ness along all the viscous boundaries can be adjusted at 
once. The hydraulically smooth k£ value of 1.0 is the 
default value. In specifying the surface roughness for 
either input, values greater than zero are treated as k£ and 
values less than zero are treated as kR (average sand-grain 
surface roughness nondimensionalized by the reference 
length). According to reference 23, the effective sand- 
grain surface roughness is classified as follows: 

k£ < 5 Hydraulically smooth 

5 < k£ < 70 

70 < k£ < 400 Transitional 


The NPARC code computes the k£ value specified for the 
individual boundary through the auxiliary pressure variable 
and compares it with the smooth-wall value given by 
SRDEF. It then uses the larger of these values to compute 
co at the boundary. Failure to specify the auxiliary pressure 
will result in a kR value of zero (default), and computations 
will be performed using the default surface roughness 
given by SRDEF. Because boundary conditions 62 and 68 
use the auxiliary pressure variable for other purposes, the 
surface roughness given by SRDEF is used. No calculations 
have been performed using any surface roughness value 
other than the default value (kR = 1.0), and users are 
advised to do the same until such time as this option can 
be more fully investigated. 

All other boundary conditions remain unchanged. As 
with the k-e model, the turbulent quantities are extrapolated 
from the interior of the flow field along slip walls and at 
free boundaries. Values of turbulence intensity I and 
turbulent viscosity p t may be specified for free inflow 
boundaries as described by Georgiadis, Chitsomboon, 
and Zhu. 20 Using a fixed inflow will cause the turbulent 
quantities along that boundary to remain unchanged from 
those in the restart file. 

Stability Considerations 


The limiters used by NPARC to increase convergence 
and stability by capping the values of the turbulent quantities 
at both the high and low extremes were modified to 
accommodate the k-co model. These limiters are used on 
the interior of the flow field, not along boundaries. The 
current implementation is similar to that used by the k-e 
model, but because co approaches infinity near a smooth 
viscous boundary, numerical stability problems may arise 
if it is limited in this region. In regions where the turbulent 
kinetic energy is found to be very small, k is set to a 
minimum value, but co is not changed. If co becomes too 
small, both k and co are set to minimum values. Should the 
turbulent kinetic energy exceed the maximum value given 
by k max = 0.10xUref(k-e)> it is set to the maximum value 
and co is computed using 

co = Maximum(co - Calculated; a* P^max^unax) 

where |it,max is specified through the user input TMUM AX. 
No upper limit check is made on co. Users who experience 
convergence difficulties or erroneous results should check 
to make sure the upper limits for the turbulent kinetic 
energy (adjusted through the variable UREFKE) and the 
turbulent viscosity (adjusted through TMUMAX) are set 
sufficiently high. The output file (26) lists the iteration 
number and the number of points where the turbulent 
quantities are limited by k^n, k^, and 00^. 
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We found that near regions of rapid acceleration or 
deceleration, such as near leading edges and behind shock 
waves, the turbulent viscosity can become very large. The 
variable option in the k-e model has been modified for 
use in the k-co model and has helped to alleviate this 
problem. The modification involves scaling the turbulent 
viscosity based on the ratio of production to dissipation R 
as: 


0.10738 (0.64286 + 0.19607R) pk 
M ' t ~ 0.09 [l + 0.357(R-l)f ® 

Details of the formulation may be found in reference 24. 
Even though the variable Cp. option may be turned on or 
off in the k-e model by adjusting the value of ICMU, ithas 
been included in the k-co model regardless. Results for the 
cases discussed next were not adversely affected by this 
modification. In fact, some results showed a slight 
improvement. 


Example 1 


Consider the two-dimensional flow over a smooth flat 
plate. For a 1 1 1x81 grid the NPARC2D formatted input 
file is given as 


&TURBIN 

IMUTUR=2, 

NTURB=5200, 

&END 

&BLOCK 

INVISC(1)=1, 

INVISC(2)=1, 

&END 


IFMAX=1, 

IMUTR2=7, 


LAMIN(1)=0, 
LAMIN(2) =1, 


NBCSEG=5, 


1 

15 

1 

1 

50 

1 



16 

111 

1 

1 

60 

1 



1 

111 

81 

81 

50 -1 



111 

111 

2 

80 

0 -1 

0.7143 

1.0000 

1 

1 

2 

80 

0 

1 

0.7345 

1.0080 


In this example, the k-co model will be initialized from 
the Baldwin-Lomax model at iteration 5200 (assuming 
the current iteration NC is less than 5200). The first 1 5 grid 
points are treated as a slip wall to allow for a uniform 
profile at the leading edge of the flat plate. The inflow and 
outflow boundaries are both free boundaries. The free 
stream is modeled using a slip boundary that is far from the 
flat plate so as not to disturb the boundary layer. The flat 
plate is smooth, and the default surface roughness is used. 
Since no compressibility corrections have been specified, 
the default setting (Sarkar’s correction to the k-equation) 
will be used. 


Example 2 


Next, consider a two-dimensional channel flow to be 
computed on a 101x5 1 grid. The formatted inputs for this 
flow are given as 


&INPUTS 

NC=5000, 

&END 

&TURBIN 

IMUTUR=4, 

COMPK=1.5, 

SRDEF=1.0, 

&END 

&BLOCK 

INVISC(1)=1, 

INVISC(2)=1, 

&END 


IMUTR2=7, 

COMPW=1.5, 


LAMIN(1)=0, 

LAMIN(2)=1, 


NTURB=5200, 

CMT0=0.25, 


NBCSEG=5, 


1 l 2 50 0 1 0.7000 1.0000 

101 101 2 50 0 -1 0.6000 0.9800 

1 31 1 1 60 1 

32 101 1 1 60 1 50.0000 

1 101 51 51 61 -1 -0.0001 0.9900 


In this case, the k-co model will be initialized from the 
k-e solution at iteration 5200. To do so, IMUTUR is set to 
the k-e model and the current iteration is reset to a value 
lower than NTURB (on subsequent runs, NC should be set 
to -1). If NC is not set below NTURB for the initial run, 
NPARC will read in k, e, and p t believing them to be k, co, 
and p t . However, the values of £ and co can differ by 
several orders of magnitude, and this may result in poor 
convergence or erroneous results or both. The com- 
pressibility correction of Wilcox has been chosen to 
override the Sarkar default settings. The minimum 
allowable surface roughness SRDEFhas been set to a kR 
value of 1 (which is what the default value would have 
been had this line been omitted). The lower boundary has 
been specified as a no-slip adiabatic surface. Since no 
surface roughness has been specified for the first 31 grid 
points, the default roughness given by SRDEF will be 
used. Along the remainder of the lower wall, a k£ value of 
50.0 will be used. The upper boundary is a no-slip iso- 
thermal wall. Along this boundary, a surface roughness 
kR (nondimensionalized by the reference length) of 0.000 1 
will be used. If the kR value corresponding to this kR at any 
point along this boundary becomes less than the smooth- 
wall surface roughness given by SRDEF (kR = 1.0), the 
value of SRDEF will be used instead. 
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Results 

Smooth Flat Plate 

Incompressible flow over a smooth flat plate was used 
as an initial test case for the k-co model. A 1 11x81 grid was 
used to model a Mach 0.2 flow. The first 15 grid points 
were treated as a slip wall in order for the flow to reach the 
leading edge of the flat plate with a uniform profile. The 
grid was sufficiently packed in the stream wise direction to 
resolve the flow gradients at the leading edge of the flat 
plate and normal to the surface to resolve the boundary 
layer. The k-co solution was initiated from the Baldwin- 
Lomax algebraic model by assuming that the production 
in turbulent kinetic energy was equal to the dissipation 
rate. The solution was determined to be converged when 
fluctuations in the mean flow and turbulent quantities at 
the last streamwise coordinate ceased. This corresponded 
to a drop in the L2 residual error of four orders of 
magnitude. 

Figure 1 shows the variation in skin friction coefficient 
Cf along the flat plate. Both the k-£ and k-co models give 
similar results and overpredict the skin friction as compared 
with the experimental data of Wieghardt. 25 Figures 2 to 
5 show the turbulent quantities and velocities at Re x 
values of lxl 0 6 , 4xl0 6 , and lxl 0 7 . From the turbulent 
kinetic energy profiles presented in figure 2, the maximum 
value predicted by the k-co model is less than that of the 
k-£ model. This was also true in the investigation conducted 
by Lang and Shih. 1 8 Figure 3 indicates close agreement in 
the turbulent viscosity values predicted by each model 
close to the wall. The large differences between the 
maximum values do not have a dramatic effect on the 
turbulent shear stress, as shown in figure 4. The small 
discontinuity in the Baldwin-Lomax results is caused by 
the turbulent viscosity, which is not smooth at the location 
where the inner and outer layers of the model meet. The 
velocity profiles presented in figure 5 compare well with 
the k-£ model and the experimental data of references 26 
and 27. Results from the three-dimensional code (fig. 6) 
were nearly identical to the two-dimensional results, except 
for a small deviation near the leading edge of the flat plate. 

At no time while running this test case was co observed 
to be limited near the viscous wall. Both k and co reached 
minimum allowable values in the free stream, as should be 
expected. The large difference in co- values of the slip wall 
and flat plate solutions did not present any numerical 
stability difficulties. 

Fraser Diffuser 

Having found that the k-co model works as well as the 
k-£ model for the flat plate, the model was then examined 
for an adverse pressure gradient test case. The conical 


diffuser of Fraser 28 represents one of the simplest of these 
types of flows. This flow enters a reducing section after 
which it passes through a straight pipe (fig. 7(a)). Here the 
subsonic flow is tripped so that it becomes fully turbulent 
before entering the 5° half-angle diverging channel. 
Experimental results indicate that the flow very nearly, 
but does not quite, reaches separation at the diffuser exit. 

A computational grid with 121 axial and 71 radial 
points was used to model this axisymmetric flow 
(fig. 7(b)). A small circular arc was used as a transition 
between sections to promote grid orthogonality. The grid 
was also clustered near the walls and the inflow to resolve 
gradients in these regions. The inflow was modeled using 
a free boundary, and the outflow was specified through the 
use of the mass-flux boundary condition. Both the k-£ and 
k-co models were initialized from a Baldwin-Lomax 
solution. The three-dimensional k-co results were initialized 
by mapping the Baldwin-Lomax axisymmetric solution 
onto a three-dimensional grid. Fifteen planes were used in 
the circumferential direction to model one quarter of the 
flow. This test was performed to validate the ability of the 
three-dimensional code to model an axisymmetric flow. 

Figure 8 compares the velocity profiles of each of the 
models at two axial locations. The skin friction along the 
diffuser wall is given in figure 9. Both the Baldwin-Lomax 
and Baldwin-Barth models predict a separation beginning 
roughly halfway through the diffuser. The k-E model 
shows little sensitivity to the adverse pressure gradient 
and predicts a very well attached flow. The k-co model, 
however, is much more accurate at predicting the near 
separation at the diffuser exit. Results from the three- 
dimensional NP ARC code were found to be nearly identical 
to the two-dimensional results. 

Sajben Diffuser 

The Sajben 22 diffuser weak- and strong-shock cases 
were selected as the next validation cases for several 
reasons: (1) They furnish an additional test for the k-co 
model in an adverse pressure gradient. (2) They are 
compressible flow cases and provide an opportunity to 
investigate the compressibility corrections of Sarkar and 
Wilcox. (3) They involve a shock-induced separation and 
are more indicative of the types of flows many NPARC 
users investigate. 

The Sajben geometry is shown in figure 10(a), and the 
81x51 grid on which both cases were computed is shown 
in figure 10(b). Both the inflow and outflow boundaries 
were specified as free boundaries. Computations on the 
weak-shock case were performed first with the k-E and 
k-co models initialized from the Baldwin-Lomax solution. 
The outflow pressure of the k-£ and k-co solutions was then 
lowered to obtain the corresponding strong-shock results. 
The three-dimensional solutions were initialized by 
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copying the corresponding two-dimensional solution onto 
consecutive planes. Comparisons of the three-dimensional 
and two dimensional results are made using the middle 
plane. 

Weak Shock .— Figure 11 shows the pressure 
distributions along the top and bottom surfaces of the 
diffuser. The results from the various models are basically 
the same with the exception of the shock location. The 
k-e model predicts the earliest shock, followed by the 
Baldwin-Lomax and k-co models. Both the k-£ and 
Baldwin-Lomax models compare well with the 
experimental data, while the k-to model predicts the shock 
location slightly further downstream. 

Figure 12 displays the velocity profiles at two locations 
downstream of the shock. From these profiles, it can be 
seen that the k-£ model predicts a larger core velocity, 
while the k-co model is more accurate in the near-wall 
regions. A comparison of the results obtained using the 
two- and three-dimensional versions of NPARC is given 
in figure 13. Both give similar results. 

Strong Shock . — The top and bottom wall pressure 
distributions for the strong-shock case are given in fig- 
ure 14. These results show that downstream of the shock, 
the Baldwin-Lomax model has difficulties converging to 
a steady-state solution for this case, and the results presented 
represent an instant in time. Examining the two-equation 
models reveals that the k-co model still predicts the shock 
location further downstream than the data suggest, but it 
does a much better j ob predicting the pressure distribution 
in the separation region. 

Velocity profiles downstream of the shock location are 
given in figures 15(a) to (c). Although none of the models 
accurately predicts the velocity profile near the upper 
wall, the k-co model provides the best agreement. It also 
gives very good results along the lower boundary. Both 
the k-e and k-co models predict roughly the same peak 
velocity, but the k-co model more accurately matches the 
core-flow profile. Results from the two- and three- 
dimensional codes were again found to be nearly identical 
(fig. 16). 

A comparison made of the different compressibility 
corrections to the k-CO model is given in figure 17. Adding 
the Sarkar correction to only the k-equation gives the best 
results. This setting was selected as the default setting. 
Using the correction of Wilcox yields nearly the same 
results as not using any compressibility correction. This 
occurs because the turbulent kinetic energy is such that the 
turbulent Mach number at all but a very few points is 
below the cut off Mach number given by CMTO. 


Conclusions and Future Work 

The k-co model was recently implemented in the NPARC 
code. Results from the validation cases indicate that this 
model is better suited to adverse pressure gradient flow 
calculations than are other turbulence models currently 
available in NPARC. A description of the model equations 
and boundary conditions was given along with a discussion 
of several other model features not completely validated 
in this study. Future work using this model will include 
further validation and an examination of the model’s 
ability to simulate transition and surface roughness. 
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Skin friction coefficient, Cf 



Figure 1 . — Local skin friction coefficient along flat plate. 
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Distance above surface, Y/8 Distance above surface, Y/8 _ Distance above surface, Y/8 





Reynolds stress, -uv/Cu*) 2 


Figure 4.— Turbulent shear stress along flat plate for Reynolds Figure 5.— Velocity profiles along flat plate for Reynolds 

number based on plate location Re x . (a) Re x = 1 xl 0^. number based on plate location Re x . (a) Re x = 1 xl 0 . 

(b) Re x = 4x1 0 6 . (c) Re x = 1 xl 0 7 . ( b ) Re x = 4x1 ° 6 - (°) Re x = 1 x1 ° 7 * 
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Skin friction coefficient, C f 



Figure 6.— Comparison of two- and three-dimensional k-to skin friction results for flat plate. 


/— Reducing section 



Straight pipe 



Figure 7.— Fraser diffuser, (a) Geometry, (b) Computational mesh. 
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Figure 8.— Fraser diffuser profiles at two axial locations x. (a) x = 1 17 mm. 
(b) x = 642 mm. 
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Figure 10. — Sajben diffuser, (a) Geometry, (b) Computational mesh. 
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Axial position, x/H 


Figure 11. — Pressure distributions for Sajben diffuser (weak shock), (a) Top surface, 
(b) Bottom surface. 
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Figure 1 2.— Sajben diffuser velocity profiles (weak shock) at two axial positions x/H. 
(a) x/H = 4.6. (b) x/H = 6.3. 
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Pressure ratio, P/P t Pressure ratio, P/P t 




Figure 14. — Pressure distributions for Sajben diffuser 
(strong shock), (a) Top surface, (b) Bottom surface. 





Figure 15. — Sajben diffuser velocity profiles (strong shock) 
at three axial positions x/H. (a) x/H = 4.6. (b) x/H = 6.3. 

(c) x/H = 7.5. 
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Vertical position, y/H Vertical position, y/H 



Figure 1 6. — Comparison of k-w results for Sajben diffuser (strong shock) at axial position x/H = 7.5. 



Figure 17.— Comparison of compressibility corrections for Sajben diffuser (strong shock) at 
axial position x/H = 6.3. 
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